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Abstract 


C" A  summary  of  the  bubble  pulsation  dynamics  of  an  underwater  explosion  is  presented 
with  the  aim  of  establishing  the  whipping  response  of  a  nearby  surface  ship.  The  effects  of 
the  initial  shock  front  are  not  considered  and  the  fluid  is  assumed  to  be  an  incompressible 
medium.  The  gas  bubble  volume  dynamics  can  then  be  related  to  far  field  fluid  accelerations 
assuming  that  a  potential  flow  condition  is  satisfied.  Given  the  effective  mass  of  the  ship 
comprising  the  added  fluid  and  bouyancy  mass  contributions,  the  forces  on  the  ship  resulting 
from  the  accelerated  flow  field  can  be  determined.  The  theory  considers  the  effects  of 
proximity  of  the  explosion  bubble  to  a  free  surface,  vertical  migration  of  the  bubble  due  to 
the  hydrostatic  imbalance,  and  an  energy  dissipative  mechanism  introduced  via  a  pseudo 
pressure  drag  coefficient.  The  theory  is  applied  in  analyses  of  the  whipping  displacements 
and  stresses  induced  in  the  hull  of  a  typical  frigate-size  warship  subjected  to  underwater 
blast  loading.  The  response  is  obtained  using  a  modal  superposition  of  the  wet  hull  modes 
of  a  detailed  equivalent  beam  model  of  the  hull.  CCn m  a  <■( a)  .  <1 


REsumE 


On  resume  la  dynandque  des  oscillations  de  la  bulle  produite  par  una 
explosion  sous-marine  en  vue  d'Etablir  las  caractEristiques  de  fouettement 
vertical  d'un  navire  voguant  non  loin  en  surface.  On  ne  tient  pas  eompte  des 
effets  du  front  de  choc  initial  et  on  suppose  que  le  fluide  est 
incompressible.  On  peut  ensuite  Etablir  un  lien  entre  la  dynamique  du  volume 
de  la  bulle  de  gaz  et  les  accelerations  du  fluide  dans  le  champ  lointain  si  on 
suppose  qu'une  condition  de  champ  de  potentiel  est  satisfaite.  Etant  donnE  la 
masse  effective  du  navire,  y  comprise  les  apports  de  la  masse  de  fluide 
ajoutee  et  de  la  masse  de  fluide  dEplacEe,  on  peut  determiner  les  forces 
exercees  sur  le  navire  par  le  champ  d'ecoulemant  accElErE.  La  thEorie  tient 
eompte  des  effete  de  la  proximite  entre  la  bulle  et  une  surface  libre,  de  la 
migration  verticals  de  la  bulle  k  cause  du  desequilibre  hydrostatique  et  d'un 
mEcanisme  de  dissipation  de  l'Energie  represent^  par  un  pseudo-coefficient  de 
trainee  de  pression.  La  thEorie  est  appliquEe  4  des  analyses  des  mouvements 
de  fouettement  vertical  et  des  contraintas  indults  dans  la  coque  d'un  bfltiment 
de  guerre  type  de  la  taille  d'une  frigate  par  une  explosion  sous-marine.  La 
rEponse  est  obtenue  4  l'aide  d'une  superposition  modale  des  modes  humides  d'un 
module  Equivalent  de  la  coque  ou  cette  demiEre  eat  reprEsentEe  en  dEtail  par 
une  poutre. 
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Notation 


a  bubble  radius 

a  time  derivative  of  bubble  radius 
A  projected  area 

Ci  drag  coefficient 

d  depth  of  charge 

ei,ci  source  and  dipole  strength  coefficients 
Eq  total  energy  of  explosion 

Ff  frequency  of  first  bubble  pulse 

g  acceleration  of  gravity 

k  non-dimensional  energy  parameter 
k\  adiabatic  constant 

K,j  keel  shock  factor 

L  length  scale  factor 

P  general  point  in  the  fluid  domain 
r  radial  distance  from  bubble  center 
t  time 

T  time  scale  factor 

Tp  period  of  first  bubble  pulse 
u  fluid  velocity 

U  kinetic  energy  of  bubble  pulse 

v  volume  of  bubble 

V  potential  energy  of  bubble 
w  weight  function 

v 


W  charge  mass 

x  non-dimensional  bubble  radius 
xo,  xm  minimum  and  maximum  bubble  radii 
z  pressure  head 

so  initial  pressure  head 
a  migration  control  coefficient 

jff  free  surface  effect  control  coefficient 

7  adiabatic  gas  constant 

S  non-dimensional  depth 

(  energy/mass  (2.051  *  10®  joules/gm  TNT  ) 
£  non-dimensional  pressure  head 

£o  non-dimensional  initial  pressure  head 

0  general  angle 

A  non-dimensional  time  derivative  of  z  (f) 
v  vertical  bubble  velocity 
p  density  of  fluid 

pg  density  of  grs 

a  non-dimensional  time  derivative  of  *  (£) 

r  non-dimensional  time 

4>  velocity  potential 

X  period  ratio  function 
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1  Introduction 


An  underwater  explosion  interacts  with  the  surrounding  fluid  in  two  different  phases 
which  are  potentially  damaging  to  a  nearby  ship  or  submarine.  The  first  phase  is  a  transient 
shock  front  which  causes  a  rapid  rise  in  the  fluid  velocity,  and  large  inertial  loading.  This 
phase  imparts  little  momentum  to  the  fluid  because  of  its  extremely  short  duration.  The 
second  phase  in  the  explosion  is  a  radial  pulsation  of  the  gas  bubble  with  a  duration  much 
longer  than  the  shock  phase.  The  pulsations  are  a  result  of  the  imbalance  of  internal 
explosion  product  pressure  and  the  fluid  hydrostatic  pressure.  Depending  on  the  depth  and 
size  of  the  charge,  several  pulsations  are  possible,  although  the  assumptions  and  theory 
presented  in  this  report  are  considered  valid  for  only  the  first  complete  pulse. 

The  bubble  pulsations  lead  to  quite  significant  pressure  impulses  on  a  nearby  ship  hull. 
As  well  as  causing  local  damage,  they  tend  to  excite  the  ship  in  heave  and  natural  vertical 
vibration  modes.  Because  the  pulsation  periods  are  often  close  to  the  two  and  three-noded 
bending  vibration  frequencies  of  a  typical  frigate-size  warship,  a  near  resonant  condition 
can  exist  with  the  potential  to  cause  large  amplitude  whipping  displacements.  These  dis¬ 
placements  in  turn  can  overstress  the  hull  and,  in  severe  cases,  lead  to  an  overall  hull  failure. 
The  factors  determining  the  whipping  response  of  a  particular  ship  to  an  underwater  blast 
load  thus  depend  on  the  characteristics  of  the  charge  as  well  as  the  vibration  characteristics 
of  the  ship. 

Explosion  bubble  dynamics  have  been  treated  by  aeveral  authors^'2’*'4),  and  the  ap¬ 
plication  of  the  theory  in  the  analysis  of  ship  whipping  has  been  discussed  extensively  by 
Hicks(6).  The  latter  work,  although  covering  many  aspects  of  the  complete  problem,  lacks 
detail  and  consistency  in  the  derivation  of  the  bubble  dynamics  and  solution  procedures 
required  in  the  application  of  the  theory.  Several  authors  have  reviewed  the  work  as  pre¬ 
sented  by  Hicks  and  others^8'7),  and  have  clarified  certain  aspects  of  the  theory;  however, 
the  derivation  and  solution  of  the  bubble  dynamical  equations  has  remained  somewhat  con¬ 
fusing.  In  this  report,  the  bubble  energy  equations  are  summarized  clearly  in  a  manner 
intended  to  facilitate  a  numerical  solution  procedure.  The  assumptions  are  explicitly  stated 
for  the  various  cases  treated,  and  a  solution  procedure  is  discussed.  The  algebraic  develop¬ 
ment  is  not  given  in  detail,  and  the  interested  reader  is  referred  to  References  2,  5  and  6 
for  further  discussion, 

The  presentation  of  the  bubble  dynamic  equations  begins  with  the  simplest  case  of  a  non¬ 
migrating  bubble  far  from  the  free  surface.  In  this  and  all  subsequent  derivations,  a  potential 
flow  regime  is  assumed  to  be  valid.  For  thU  first  case,  the  single  differential  equation 
describing  the  bubble  radius  can  be  effectively  integrated  with  Tschebycheff  polynomials. 
The  introduction  of  bubble  migration  and  free  surface  effects  result  in  a  system  of  first 
order  differential  equations  which  can  be  integrated  via  a  Runge-Kutta  algorithm.  A  drag 
term  analogous  to  that  for  a  solid  body  wake  can  be  introduced  in  the  bubble  migration 
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equation  to  account  for  energy  dissipation  during  the  bubble  expansion  and  contraction. 
This  term  has  little  theoretical  justification,  but  provides  an  adjustment  parameter  to  give 
energy  dissipation  rates  in  agreement  with  experimental  results. 

The  solution  of  the  bubble  expansion  and  migration  equations  provides  an  analytical 
description  of  the  velocity  potential.  The  far  field  fluid  velocities,  accelerations,  and  pres¬ 
sures  can  then  be  obtained  from  application  of  the  appropriate  operators  on  the  velocity 
potential  in  accordance  with  potential  flow  theory  and  general  fluid  dynamics.  The  fluid 
accelerations,  which  are  used  in  the  current  approach,  define  a  time-dependent  forcing  func¬ 
tion  per  unit  mass  at  a  point  in  the  far  field.  This  function  can  be  spatially  averaged  and 
scaled  by  the  appropriate  mass  term  to  provide  a  point  load  force  time  history  representing 
the  effective  load  on  a  hull  segment.  The  response  of  the  ship  can  then  be  determined  via 
direct  integration  methods  or,  more  commonly,  by  modal  superposition. 

Because  of  the  simplicity  and  accepted  accuracy  of  equivalent  beam  representations 
of  ships  in  the  evaluation  of  low  frequency  vibrations,  such  models  are  often  used  in  the 
prediction  of  whipping  response  to  underwater  explosions.  Typically,  a  minimum  of  twenty 
sections  are  used  to  represent  both  the  hydrodynamic  form  and  structural  characteristics  of 
the  ship.  A  recent  study  indicates  that  this  is  an  acceptable  level  of  discretization  for  the 
determination  of  the  lowest  modes  of  vibration  of  a  frigate  hull.!8)  In  the  present  study,  one 
hundred  sections  are  used  for  the  structural  representation  of  a  test  ship  typical  of  a  modern 
frigate.  Twenty  sections  are  used  for  the  hydrodynamic  description.  The  hydrodynamic 
coefficients  of  form  can  be  used  with  the  traditional  Lewis  form  methods!0)  to  obtain  the 
mode-dependent  sectional  added  fluid  masses  associated  with  hull  vibration.  Alternative 
methods  of'  obtaining  added  masses  independent  of  the  vibration  mode  are  also  available, 
and  are  applied  in  the  determination  of  the  wet  modes  of  the  hull  in  the  current  study. 
Once  the  added  mass  terms  are  known,  an  eigenvalue  analysis  of  the  beam  model  including 
this  mass  will  provide  with  sufficient  accuracy  the  first  S  or  4  vertical  bending  modes  of 
the  hull.  These  modes  in  turn  define  a  suitable  basis  for  solution  of  the  forced  equations  of 
motion  via  modal  superposition. 

Section  2  of  this  report  presents  the  derivations  of  the  bubble  equations,  and  Section 
3  discusses  the  application  of  the  theory  in  the  analysis  of  the  response  of  the  test  ship  to 
underwater  blast  loads  using  the  methods  outlined  in  the  text.  A  final  section  includes  some 
concluding  remarks  on  the  reliabiity  of  the  theory  and  general  comments  on  the  results  of 
the  analysis. 

2  Bubble  Dynamics 

A  number  of  assumptions  are  implicit  in  the  derivation  of  the  bubble  pulsation  and 
migration  equations,  and  the  prediction  of  the  response  of  nearby  bodies  to  the  bubble 


dynamics.  First,  as  discussed  above,  the  compressible  phase  of  the  explosion  is  neglected 
completely  on  the  basis  that  the  impulse  of  that  phase  is  very  small  compared  to  the 
pressure  impulses  of  the  bubble  pulsations.  Secondly,  the  fluid  is  considered  incompressible 
and  irrotational,  and  a  potential  flow  regime  is  assumed  to  represent  the  fluid  motion  caused 
by  the  expansion  and  contraction  of  the  bubble.  Thirdly,  the  bubble  shape  is  arsumed  to 
remain  spherical  throughout  the  pulsations.  This  is  an  adequate  approximation  for  the 
expansion  phase  although  large  distortions  tend  to  occur  near  the  end  of  the  contraction 
phase.  These  deviations  from  the  spherical  shape  make  this  assumption  a  poor  one  for 
the  analysis  of  subsequent  bubble  pulses.  Finally,  a  far  field  approximation  is  made  in  the 
derivations  in  this  report  whereby  the  ship  is  assumed  to  have  no  influence  on  the  bubble 
dynamics.  The  effects  of  close  proximity  and  contact  explosions  can  be  treated;^)  however, 
in  that  case  local  damage  rather  than  whipping  is  the  primary  concern. 

Several  cases  are  treated  in.  order  of  complexity  beginning  with  the  simple  case  of  a 
non-migrating  bubble  in  deep  water.  In  all  cases  an  energy  conservation  approach  is  used 
to  derive  the  pulsation  dynamics,  and  the  losses  due  to  acoustic  radiation  and  cooling  of  the 
gas  products  are  neglected.  This  last  point  is  very  important  since  it  essentially  limits  the 
accuracy  of  this  form  of  derivation  to  single  bubble  pulses.  The  energy  dissipative  mech¬ 
anisms  which  control  the  inter-bubble  phase  are  not  well  understood,  and  approximations 
are  all  that  is  possible  in  this  approach. 

2.1  Deep  bubble  approximation 

At  significant  depth  and  hence  large  ambient  pressure,  an  underwater  explosion  bubble 
is  not  affected  by  the  free  surface  reflection,  and  undergoes  little  vertical  migration  because 
of  the  low  bouyancy  forces  associated  with  a  reduced  bubble  radiuB.  The  bubble  center  can 
be  considered  as  a  stationary,  singular  source  with  a  time  dependent  strength.  The  fluid 
velocity  potential  resulting  from  the  source  becomes, 


The  time  dependent  source  strength  parameter,  ei,  can  be  determined  by  considering  the 
fluid  velocity  at  the  surface  of  the  bubble.  The  bubble  volume  time  rate  of  change  is, 


hence, 
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The  fluid  velocity  at  a  distance  r  from  the  explosion  source  is  then, 


u  = 


a2a 


(5) 


An  energy  equation  based  on  the  conservation  of  total  energy  in  the  system  can  be  formu¬ 
lated, 

E0  =  U  +  V  (6) 

where  Eq  is  the  initial  energy  of  the  explosive,  determined  from  the  knowledge  of  the  mass 
and  type  of  charge,  and  U  and  V  are  the  kinetic  and  potential  energies  of  the  explosion 
bubble.  The  kinetic  energy  associated  with  the  explosion  bubble  is  the  energy  of  the  Bur- 
rounding  fluid,  given  by  the  integral, 

U  ~\ja  P^H***)  dr  (7) 

The  potential  energy  is  the  sum  of  the  pressure  potential  of  the  explosion  void  and  the 
internal  energy  of  the  gas  products.  The  complete  energy  equation  can  be  expressed  in 
terms  of  the  bubble  radius,  a,  and  its  derivatives  as,(*) 

r  -  ,»,■  ,  ,  *1  Wl’'1 

E„  -  2,  pa  a  +  -a pa  g,„  +  (8) 


where  an  adiabatic  equation  of  state  has  been  assumed  valid  for  the  explosion  reaction. 
This  equation  can  be  written  in  the  much  simpler  form, 
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by  introducing  the  non-dimensional  factors, 


with, 
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(11) 

(12) 

(13) 

(14) 


Based  on  experimental  evidence/10)  it  is  necessary  that  7  —  f  if  equation  (9)  is  to  describe 
the  behaviour  of  the  bubble.  For  the  theoretical  and  experimental  bubble  half  periods,  T, 
to  agree,  the  following  relation  is  specified, 

E0  =  tW  (15) 

with  e  as  the  energy /unit  mass  given  as  approximately  2.051  *  10®  joules/gram  for  TNT. 
Comparisons  of  peak  pressures  indicate  that, 

k  *O.0743(zo)«  (16) 

for  units  of  meters,  kilograms  and  seconds.  Hence  the  specification  of  the  charge  type,  mass 
and  depth  allows  the  calculation  of  the  time  and  length  scale  factors,  and  it  remains  to  find 
the  solution  x(r),  given  by  equation  (9).  This  equation  can  be  written  as, 


/x-3  -  fcE-37  -  1 


Integrating, 

_  /•*(*■)  dx 

J x(r=0)  \/x~ 3  -  k»~l'1  -  1 

This  integral  can  be  evaluated  very  effectively  using  TschebychefF  polynomials/11)  where 
the  integral  I  defined  as, 

r_  fk _ fjy) _  j., 


-  a  0- 


can  be  integrated  as, 


/=!>/(*) +  *n 


where  fin  is  a  small  remainder  term.  For  the  integrals  of  interest  in  the  present  study, 
near  exact  integration  is  obtained  using  an  eight  term  expansion  (n  =  8).  The  function  is 
evaluated  at, 

*  =  <*±0  +  ^*  (21) 


(2*  -  1) 
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and  the  weight  function  is  in  this  case  a  constant, 


Using  this  approach,  the  bubble  radius  time  function  is  rewritten  as, 


fXm  V(*  ~  a)(&  ~ 
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(24) 


It  is  evident  that  the  evaluation  of  this  integral  provides  the  time  at  which  the  radius 
assumes  a  particular  value,  as  opposed  to  a  more  common  time  integration  process.  The 
integrand  causes  problems  when  the  denominator  function  goes  to  zero.  This  occurs  at  the 
minimum  and  maximum  radii  of  the  bubble,  *o  and  xm  respectively,  and  these  values  must 
become  the  integration  limits  to  avoid  singularities  in  the  integral.  These  points  can  be 
determined  from  the  roots  of  the  denominator  function,  which  is  equivalent  to  setting  x  =  0 
in  equation  (9).  In  fact,  the  roots  of  this  expression  are  quite  sensitive  to  the  value  of  Jfc, 
which  reflects  the  dependence  of  the  bubble  dynamics  on  the  energy  of  the  bubble  and  the 
ambient  pressure.  In  comparisons  with  small  scale  experiments  from  which  the  theory  is 
calibrated/2'6)  the  magnitude  of  k  remains  small,  (of  order  0.1  •  0.2).  In  Buch  cases,  this 
approach  works  well  even  for  multiple  pulses;  however,  for  realistic  depths  of  discharge  in 
the  ship  whipping  context,  the  k  value  is  initially  quite  low,  but  increases  rapidly  beyond  a 
‘small’  value  when  energy  losses  are  introduced  in  attempts  to  apply  the  theory  to  multiple 
bubble  pulses.  In  many  cues,  no  positive  roots  of  the  function  will  exist  for  other  than  the 
first  pulse,  and  the  method  cannot  predict  multiple  pulse  behaviour.  If  reasonably  small 
energy  attenuations  are  assumed,  multiple  pulses  can  be  predicted  by  the  theory. 

The  description  of  energy  attenuation,  which  is  critical  in  the  application  of  thiB  ap¬ 
proach  for  prediction  of  multiple  bubble  pulses,  is  not  well  defined.  Classical  literature 
suggests  losses  of  approximately  70  percent  of  the  initial  energy  for  the  second  pulse,  and 
a  further  50  percent  reduction  for  the  third  pulse.  Implementation  of  these  values  in  the 
above  derivation  leads  to  only  one  bubble  pulse  for  realistic  charge  depths.  Other  energy 
attenuations  might  be  used,  since  results  of  small  scale  experiments  indicated  reasonably 
constant  reductions  of  roughly  60  percent  for  each  pulse.  These  latter  reductions  are  used 
to  obtain  the  results  presented  in  Section  3. 

2.2  Effects  of  migration 

The  vertical  migration  of  the  gas  bubble  can  be  included  very  simply  in  the  definition 
of  the  velocity  potential  function.  In  this  case,  Figure  2,  the  potential  function  is  given  as, 

<f>  =  —  -f^|cosfl  (25) 

r 

The  strength  coefficients  can  be  evaluated  on  the  bubble  surface  to  give, 

ei  =  -j-  =  a2a  (26) 
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where  v  is  the  mean  vertical  bubble  velocity.  Applying  the  conservation  of  energy  procedure 
as  in  the  first  case,  the  energy  equation  becomes, 
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The  first  term  corresponds  to  the  kinetic  energy  of  the  fluid,  the  second  term  the  kinetic 
energy  of  migration,  the  third  the  gravitational  potential,  and  the  fourth  term  the  pressure 
potential  of  the  explosion  products.  This  expression  can  be  non-dimensionalized  to  give  the 
bubble  equation  of  motion  in  the  form, 


^  +  x^+x^_. 
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where  fo  is  the  initial  non-dimensionalized  total  head  at  discharge  depth,  and  f  is  the  non- 
dimensionalized  head  as  a  function  of  time.  An  application  of  Lagrange’s  equation  to  the 
energy  expressions, 

|(^-V)  =  A(E,-V)  _  (30) 

gives  the  equation  of  motion  governing  the  migration  of  the  bubble  as, 
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Differentiating,  we  can  form  a  system  of  first  order  equation  suitable  for  integration  with  a 
standard  Runge-Kutta  integration  algorithm, 


er  =  x 
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The  initial  conditions  for  the  solution  of  this  system  are  taken  as, 
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a  —  0 
A  =  0 
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where  xo  is  taken  as  the  minimum  root  of  the  radius  function  discussed  in  the  case  of 
the  non-migrating  bubble.  The  solution  of  these  equations  defines  the  velocity  potential 
function,  which  in  turn  will  provide  the  far  field  fluid  pressures  and  accelerations.  The 
derivation  of  the  fluid  accelerations  is  discussed  in  the  final  case  of  this  section. 


2.3  Free  surface  effects 


Neglecting  the  existence  of  gravity  waves,  the  free  surface  is  defined  as  a  contour  of  zero 
potential.  This  can  be  satisfied  by  defining  the  surface  as  a  plane  of  anti-symmetry,  and 
placing  a  negative  image  bubble  above  the  plane.  The  effect  of  this  image  is  a  cancelation 
of  the  horizontal  pressure  gradient  and  accelerations  at  the  free  surface,  and  a  doubling  of 
the  vertical  components.  This  can  be  done  very  simply  as  the  last  stage  of  the  solution 
procedure.  Such  an  approach  implies  that  the  free  surface  has  no  effect  on  the  bubble 
dynamics,  which  is  true  for  reasonably  deep  charges;  however,  where  the  bubble  pulse  is 
within  several  radii  of  the  free  surface,  the  presence  of  the  surface  reflection  does  influence 
the  motion  of  the  bubble,  and  the  dynamical  equations  must  be  modified.  In  this  case,  the 
velocity  potential  is  reformulated  as, 

4>  =  —  +■  %  cos  0  -  —  +  ~  cos  0i  (40) 
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Enforcing  the  bubble  surface  velocities,  the  strength  coefficients  are  found  to  be, 


ei  s=  a? a  (41) 


where  d  is  the  depth  of  the  bubble  center  below  the  free  surface,  Figure  3.  In  fact,  typ¬ 
ical  of  imaging  methods,  an  infinite  number  of  images  are  required  to  balance  the  higher 
order  terms  introduced  into  the  expression  for  the  velocity  potential,  and  the  use  of  only  a 
single  image  is  an  approximation.  Neglecting  terms  of  order  higher  than  (J)J ,  the  energy 
conservation  equation  becomes, 
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Non-dimensionalizing,  we  obtain, 
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where  6  =  Again  using  the  Lagrange  equation,  the  bubble  migration  equation  becomes, 

r3*4ia  .  *5i3  .  *31  (4g) 
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Taking  the  derivatives,  a  system  of  first  order  equations  can  be  formed  similar  to  those  of 
equations  (32)>(35),  and  the  same  initial  conditions  can  be  used. 

2.4  Bubble  drag  forces 

Experimental  results  have  shown  that  the  migration  velocities  predicted  by  the  theory 
outlined  above  are  too  large.  To  control  these  velocities,  a  pseudo  drag  coefficient  is  intro¬ 
duced  and  the  energy  dissipation  of  that  force  is  included  in  the  overall  energy  conservation 
equation.  This  drag  term  is  not  particularity  realistic,  since  drag  coefficients  reflect  the 
interaction  of  a  solid  body  and  the  surrounding  fluid,  whereas  the  explosion  bubble  is  a 
surface  of  constant  pressure.  However,  the  controllable  energy  dissipation  of  the  drag  force 
allows  a  better  correlation  of  theoretical  and  experimental  results.  The  drag  force  is  defined 
in  the  conventional  manner, 
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where  A  is  the  projected  area  of  the  bubble,  in  this  case  a  circle  of  radius  a.  Although  the 
projected  area  of  the  bubble  is  largest  at  its  maximum  radius,  the  largest  drag  forces  and 
energy  dissipation  occur  near  the  bubble  minimum  radius.  This  is  a  result  of  the  much 
higher  migration  velocities  which  are  imparted  to  the  bubble  due  to  the  momentum  in  the 
fluid  acting  on  a  body  with  a  much  decreased  added  fluid  mass.  Inclusion  of  the  drag  force 
in  the  energy  conservation  equation  leads  to  the  bubble  expansion  equation  as, 


and  the  migration  equation, 
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Performing  the  differentiation,  we  obtain  the  system  of  equations, 
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The  a  and  /3  are  convenient  programming  controls  which  can  be  set  to  a  value  of  one  if 
the  effects  of  migration  and  the  free  surface  respectively  are  to  be  included.  The  equations 
can  then  be  programmed  in  one  subroutine  for  use  in  all  cases  in  which  the  initial  conditions 
of  equations  (36)-{39)  are  utilized. 

A  drag  coefficient  can  be  determined  by  matching  the  energy  dissipation  and  migration 
predictions  to  experimental  results,  and  a  value  of  approximately  2.0  -  2.5  appears  to- give 
reasonable  agreement.  This  is  well  in  excess  of  normal  drag  coefficients  for  solid  bodies, 
hence  either  the  assumption  of  spherical  shape  is  not  realistic,  or  the  drag  theory  does  not 
apply  to  moving  pressure  surfaces.  Regardless  of  its  theoretical  basis,  the  drag  mechanism 
can  be  used  for  correlation  purposes. 

The  strength  coefficients  of  the  velocity  potential  of  equation  (40)  can  be  expressed  in 
terms  of  the  scale  factors  and  non-dimensional  functions  as, 

•l  *  tjr*  (63) 
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The  time  derivatives  of  these  coefficients  required  in  the  derivation  of  the  fluid  accelerations 
are  then  given  as, 

e'l  =  (2®crs  +  x“b)  (55) 
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2.5  Fluid  accelerations 
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Once  the  fluid  velocity  function  is  defined,  the  fluid  velocities  and  accelerations  in  the 
far  field  can  be  calculated.  For  potential  flow,  the  velocity  is  the  negative  gradient  of  the 
potential, 

u„  =  -  V  (67) 

In  cartesian  coordinates,  Figure  2,  this  gives, 
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From  the  general  theory  of  fluid  dynamics,  the  fluid  accelerations  are  given  as, 
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The  two  terms  of  this  expression  can  be  evaluated  to  give, 
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These  expressions  can  be  substituted  into  equation  (59)  to  obtain  the  far  field  fluid  ac¬ 
celerations.  Once  the  accelerations  are  known,  the  forces  can  be  calculated  by  averaging 
the  accelerations  over  segments  of  the  ship  hull,  and  scaling  by  the  effective  mass  of  the 
section.  The  mass  term  will  include  the  bouyancy  mass  of  the  hull,  and  the  added  fluid 
mass  associated  with  vertical  vibration  of  the  hull.  These  terms  can  be  found  by  various 
methods,  and  are  assumed  to  exist  in  the  formulation  given  in  this  report. 

2.6  Solution  outline 

Based  on  the  previous  development,  the  generation  of  a  load  time  history  for  an  explosion 
can  be  summarised  in  the  following  steps: 

1.  Establish  initial  conditions  and  scale  factors  of  the  bubble  pulse  from  the  specified 
depth,  charge  mass  and  energy  conversion  factor  (Equations  (9)-(14)). 

2.  Begin  integration  of  the  relevant  equations  from  the  initial  conditions.  The  integration 
time  step  should  be  some  appropriately  small  fraction  of  the  overall  period  of  the  pulse. 

3.  At  each  time  step,  calculate  ei, «2, «i» «2»  the  fluid  acceleration  u„  and  the  force  miiv 
(Equations  (53)-(5G)  and  (59)). 

4.  Finish  integration  at  a  predetermined  cut-off  time.  This  can  be  estimated  from  the 
time  scale  factor.  If  migration  is  included,  the  integration  is  terminated  when  the 
bubble  surface  contacts  the  free  surface. 
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5.  If  multiple  pulses  are  possible,  apply  the  energy  dissipation  criterion  at  the  bubble 
minimum  radii,  and  repeat  steps  1-4  above. 

Tbe  above  approach  has  been  adopted  in  the  coding  of  a  short,  interactive  computer 
program,  UBLOAD.  This  program  accesses  the  equivalent  beam  properties  data  file  created 
by  SCRAP(1J)  to  obtain  the  sectional  fluid  masses.  The  sectional  masses  are  used  in  con¬ 
junction  with  the  hull  form  coefficient  data  file  to  calculate  the  effective  mass  of  the  ship. 
Utilizing  user-specified  information  on  the  charge  depth,  standoff,  and  mass,  the  program 
calculates  the  time  history  of  the  bubble  motion  with  either  the  Tchebycheff  polynomial 
(for  the  deep  bubble  approximation  only)  or  Runge-Kutta  integration  method.  A  sample 
program  run  with  explanatory  comments  is  included  in  the  Appendix.  The  far  field  fluid 
accelerations  and  force  time  history  are  calculated  for  all  nodes  available  in  the  structural 
representation  of  the  ship.  This  data  is  stored  in  a  format  suitable  for  direct  use  in  the 
finite  element  program  VAST,(13) 

3  Predictions  from  the  theory 

The  characteristics  of  two  underwater  explosions  are  presented  to  illustrate  the  results 
of  the  methods  discussed  in  Section  2.  The  charge  characteristics  have  been  taken  such  that 
comparisons  with  results  available  in  the  literature^’6)  can  be  made.  Very  close  agreement 
with  the  available  literature  has  been  obtained  for  these  bubble  characteristic  time  histories. 

The  first  blast  corresponds  to  a  small,  shal.’  jw  charge  with  a  mass  of  45.5  Kg.  of  TNT 
at  a  depth  of  7.6  m.  Figure  4  presents  the  first  pulse  radius  time  history  for  the  various 
assumptions  available  in  the  method.  In  this  case,  the  non-migrating  results  are  derived 
from  the  Tchebycheff  method,  which  appears  to  give  a  pulse  period  which  is  too  low  in 
comparison  to  that  obtained  using  the  Runge-Kutta  methods,  Figure  5.  The  discrepancy 
appears  to  be  a  result  of  the  time  step  and  symmetrizing  algorithm  used  in  the  Tchebycheff 
approach.  Although  that  method  could  be  improved,  it  is  simpler  to  use  the  Runge-Kutta 
method  for  all  cases,  and  this  has  been  done  in  all  further  tests. 

The  bubble  radius,  depth,  velocity,  and  pressure  radius  time  histories  for  the  shallow 
charge  are  presented  in  Figures  5,6,7  and  8  respectively.  The  results  for  the  migrating  bubble 
stop  when  the  bubble  surface  reaches  the  fluid  surface,  hence  these  are  more  realistic  than 
those  from  tho  non-migrating  assumption,  which  predicts  multiple  pulsations.  It  is  evident 
that  the  maximum  bubble  radius  is  not  sensitive  to  the  various  assumptions;  however,  the 
bubble  period  changes  significantly.  The  effect  of  the  free  surface  is  to  decrease  the  bubble 
period.  Figure  7  indicates  the  phenomenon  associated  with  the  fluid  momentum  lag  which 
causes  large  peaks  in  the  bubble  velocity  time  history  at  the  time  of  the  minimum  radius. 
These  peaks  are  largest  when  the  effects  of  the  free  surface  are  included. 
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The  pressure  radius  and  maximum  fluid  acceleration  time  histories  of  Figures  8  and  9 
indicate  that  the  neglect  of  bubble  migration  will  lead  to  the  prediction  of  erroneously  large 
pressure  and  acceleration  peaks  and  multiple  pulses.  The  error  caused  by  the  neglect  of 
migration  is  not  as  serious  as  would  appear  from  consideration  of  only  maximum  values, 
Bince  the  response  of  a  body  is  dependent  primarily  on  the  impulse  of  the  forcing  function 
rather  than  simply  the  magnitude  of  the  quite  transient  peaks.  The  variation  in  the  impulses 
predicted  for  the  three  assumptions  is  in  fact  fairly  low;  however,  the  prediction  of  multiple 
pulses  in  this  case  is  not  an  acceptable  result. 

The  second  blast  modelled  with  the  theory  is  a  large,  deep  charge  of  22?  Kg,  at  a  depth 
of  45.7  m.  The  time  history  characteristics  of  this  charge  are  presented  in  Figures  10  - 
14:  The  theory  predicts  two  complete  pulses  before  the  radius  function  becomes  singular. 
The  effect  of  the  energy  dissipation  factor  between  the  pulses  is  evident  in  the  lower  bubble 
maximum  radius  of  the  second  pulse.  In  this  case,  neither  the  bubble  radius  or  period 
are  particularily  sensitive  to  the  assumptions  of  migration,  indicating  that  the  deep  bubble 
approximation  would  be  reasonably  accurate  for  this  range  of  depth. 

3.1  Bubble  pulse  frequency  specification 

Based  on  a  combination  of  physical  intuition  and  previously  obtained  results,  a  worst 
case  whipping  scenario  for  a  ship  subjected  to  a  blast  load  will  result  when  the  following 
conditions  are  met: 

1.  The  charge  is  directly  under  the  ship  (no  standoff). 

2.  The  charge  is  approximately  amidships. 

3.  The  bubble  pulse  frequency  and  first  natural  hull  bending  frequency  are  matched. 

4.  There  are  multiple  pulses. 

The  first  two  criteria  can  easily  be  satisfied  by  defining  the  appropriate  charge  location.  It  is 
not  as  evident  how  the  latter  two  conditions  can  be  fulfilled;  however,  several  relationships 
can  be  exploited  to  obtain  the  charge  characteristics  which  result  in  a  specific  bubble  pulse 
period.  From  equations  (12)  ,(13)  and  (15) ,  the  time  scale  parameter  T  is  found  to  be  related 
to  the  charge  mass  and  pressure  head  at  the  charge  depth  by  the  relation, 
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for  the  meter,  kilogram,  second  unit  system,  where, 
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The  keel  shock  factor^14)  in  this  system  is  defined 
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Based  on  a  series  of  numerical  tests,  the  time  scalo  parameter  can  be  related  in  an  approxi¬ 
mate  manner  to  the  first  bubble  pulse  period  Tp,  as  shown  in  Figure  15.  These  curves  have 
been  generated  from  repeated  analyses  with  migration  and  free  surface  effects  included.  The 
relationship  is  evidently  a  function  of  depth  but  is  not  strongly  dependent  on  the  charge 
mass  in  the  range  of  variables  appropriate  in  the  whipping  problem.  This  function  is  de¬ 
noted  x(*)>  has  been  approximated  with  three  linear  segments  for  the  purposes  of  this 
study.  We  can  then  write, 
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For  a  specified  keel  shock  factor,  the  bubble  pulse  frequency  then  becomes, 
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This  equation  has  been  used  to  generate  the  curves  shown  in  Figure  16.  The  specification 
of  a  frequency  then  allows  a  charge  mass  to  be  chosen  such  that  a  specific  keel  shock  factor 
will  result.  The  charge  depth  for  that  keel  shock  factor  can  be  obtained  from  Figure  17. 
The  keel  shock  factor  in  this  system  simply  provides  a  convenient  relationship  between  the 
charge  mass,  depth,  and  explosion  strength  relative  to  the  hull. 

It  is  evident  from  Figure  16  that  there  are  numerous  combinations  of  charge  depth  and 
mass  which  will  give  a  specific  bubble  pulse  frequency.  If  matched  to  the  lower  frequencies 
representative  of  the  primary  bending  modes  of  a  frigate-sized  hull,  these  combinations 
will  result  in  high  keel  shock  factors.  In  such  cases,  local  damage  as  well  as  whipping  will 
no  doubt  be  important.  Lower  keel  shock  factors  can  be  obtained  if  the  second  vertical 
bending  and  bubble  pulse  frequencies  are  matched.  The  selection  of  a  particular  charge 
mass  and  depth  from  Figures  16  and  17  does  not  necessarily  provide  multiple  pulses.  This 
is  difficult  to  predict  a  priori  because  of  the  effects  of  migrntion  on  the  bubble  depth  and 
period,  and  a  more  comprehensive  numerical  investigation  would  be  required  to  provide 
some  simple  moans  of  predicting  the  occurrence  of  multiple  pulses.  Since  the  theory  is 
inherently  less  accurate  for  eubsequent  pulses  scause  of  the  unrepresentative  modelling  of 
the  energy  dissipation  and  bubble  distortion  meei,  anisms,  such  an  investigation  is  probably 
not  justified. 


3.2  Application  of  the  Theory 

It  is  informative  and  indeed  the  end  objective  to  apply  the  whipping  theory  in  the 
prediction  of  the  response  of  an  actual  ship.  To  this  end,  a  very  detailed  equivalent  beam 
model  of  a  typical  frigate-size  warship  has  been  created,  from  which  the  wet  hull  modes  could 
be  calculated.  Included  in  the  prediction  of  these  natural  frequencies  were  the  effects  of 
shear  area  coefficients,  deckhouse  stiffness  contribution,  and  consideration  of  the  structurally 
ineffective  material  in  the  upper  decks.  The  added  fluid  mass  associated  with  the  hull 
vibration  modes  was  calculated  using  the  line  dipole  method  presented  by  Hicks. The 
first  few  natural  frequencies  obtained  from  this  model  are  presented  in  Table  1. 


TABLE  li  Natural  frequencies  of  the  test  ship  (Hz.) 


Mode  1" 

Mode  2* 

Mode  3" 

Mode  4° 

Mode  5‘ 

Me  f,  8* 

1.68 

2,48 

3.94 

4.62 

5.14 

6.25 

"Vertical  bending  mode 
6  Horizontal  bending  mode 
6  Coupled  bending  mode 


Using  Figures  16  and  17,  a  charge  mass  and  depth  were  selected  to  provide  a  first  bubble 
frequency  of  1.68  Hz  For  the  purposes  of  this  example,  a  keel  shock  factor  of  approximately 
1.2  was  used.  The  design  charts  indicate  that  a  mass  of  approximately  167  Kg  and  depth 
of  24  m  would  result  in  the  appropriate  bubble  frequency.  The  frequency  from  the  inte¬ 
grated  time  history  was  1.66  Hz,  indicating  that  the  design  method  can  give  quite  adequate 
results.  The  charge  was  placed  directly  below  the  keel  at  a  point  70  m  from  the  forward 
perpendicular,  a  position  selected  on  the  basis  of  a  brief  study  of  the  charge  position  versus 
hull  response  for  a  constant  depth  and  mass  chargeJ”)  For  the  test  ship,  this  charge  position 
resulted  in  the  largest  hull  displacements. 

The  hull  response  was  calculated  using  modal  superposition  of  the  first  three  vertical 
hull  modes.  The  acceleration  time  history  responses  of  the  bow  and  stern  of  the  ship  are 
presented  in  Figure  18.  The  contribution  of  higher  modes  is  evident,  but  the  overall  response 
is  controlled  by  the  fundamental  mode  with  period  0.58  seconds.  A  peak  acceleration  of 
slightly  greater  than  1.0  g  is  predicted  at  the  bow,  and  a  peak  of  approximately  0.75  g  at 
the  stem.  This  is  consistent  with  the  primary  bending  mode  shape  of  the  ship,  which  has 
larger  displacements  in  the  bow  than  the  stern.1 !*)  The  higher  modal  effects  are  not  evident 
in  the  displacement  response,  Figure  Id,  which  is  completely  dominated  by  the  primary 
bending  mode.  In  an  explosion  not  tuned  to  the  primary  bending  mode,  contributions  from 
the  second  and  third  modes  have  been  noted,  although  they  are  quite  nominal.  A  maximum 
displacement  of  approximately  73  mm  occurs  at  the  bow  due  to  this  explosion. 

The  maximum  stress  induced  in  the  hull  from  the  whipping  motions  is  predicted  to 
be  70-80  MPa  in  the  midship  area.  This  stress  level  is  well  below  the  midship  stresses 


expected  from  the  design  wave  loading^8)  and  the  response  in  this  case  is  completely  elastic. 
The  maximum  pressures  at  the  keel  are  not  expected  to  cause  any  permanent  local  plate 
deformation  for  this  load  case. 

4  Conclusion 

The  spherical  explosion  bubble  theory  is  fairly  widely  used  to  predict  the  whipping 
response  of  ship  hulls  to  underwater  blast  loads,  The  theory  as  presented  is  relatively 
simple,  if  not  completely  physically  realistic,  The  inclusion  of  the  effects  of  the  free  surface 
and  bubble  migration  appear  to  give  predictable  behaviour  in  the  bubble  characteristic  time 
histories,  and  simple  numerical  techniques  can  be  UBed  to  calculate  those  time  histories. 
The  major  drawback  in  the  theory  is,  from  the  point  of  view  of  this  author,  its  inability  to 
accurately  model  the  multiple  pulse  phenomenon;  however,  it  is  conservative  to  adopt  the 
approach  taken  in  this  study,  in  which  the  charge  masses,  and  hence  bubble  energies,  are 
scaled  from  pulse  to  pulse.  Unfortunately,  the  value  of  that  scale  factor  is  not  well  defined. 

Regardless  of  the  specific  limitations  of  the  current  approach,  the  whipping  response 
analysis  is  a  global  response  problem,  and  the  magnitude  of  hull  displacements  or  accel¬ 
erations  will  be  little  affected  by  minor  changes  in  bubble  characteristics.  The  methods 
adopted  in  the  UBLOAD  program  are  stable,  efficient,  and  give  results  consistent  with 
intuition. 

The  design  charts  developed  to  allow  tuning  of  the  bubble  pulse  and  hull  natural  fre¬ 
quency  can  bo  used  effectively  to  generate  a  worst  case  loading  scenario  for  a  given  keel 
shock  factor.  The  application  of  the  method  in  the  prediction  of  the  whipping  response  of  a 
ship  is  straightforward,  given  an  equivalent  beam  representation  of  the  hull,  and  the  hydro¬ 
dynamic  characteristics  of  the  hull  form.  The  study  of  the  frigate  response  to  underwater 
blast  loading  has  indicated  the  form  and  magnitude  of  the  response  to  be  expected  in  such 
a  ship,  and  comparisons  of  the  predicted  response  to  the  results  of  any  large  scale  shock 
trials  of  this  ship  would  be  valuable. 


Figure  1:  Geometry  of  explosion  bubble  and  far  field 


Figure  2:  Geometry  of  migrating  explosion  bubble  and  far  field  velocity 
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Figure  3:  Image  and  Bource  bubble  geometry 
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Figure  4:  Bubble  radius  time  history  for  the  first  pulse  for  a  charge  of  45.5  Kg.  of  TNT  at 
a  depth  of  7.6m. 
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Figure  6:  Bubble  depth  time  history  for  a  charge  of  45.6  Kg.  of  TNT  at  a  depth  of  7.6m. 
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Figure  10:  Bubble  radius  time  history  for  a  charge  of  227  Kg.  of  TNT  at  a  depth  of  45  m. 
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Figure  13:  Pressure-radius  time  history  for  a  charge  of  227  Kg.  of  TNT  at  a  depth  of  45  m 


28 


[ME  SCALE/ACTUAL  PERIOD  T/TP 


0,80 

0.78 

0,76 

0.74 

0.72 

0.70 

0.68 

0.66 

0  10  20  30  40  50 

DEPTH  Cm) 

Figure  15:  Period  ratio  as  a  function  of  depth  for  various  charge  masses. 
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Figure  18:  Acceleration  response  of  the  bow  and  stern  of  the  test  ship  to  a  tuned  bubble 
pulse  loading  (keel  shock  factor  1.19). 
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Figure  10:  Displacement  response  of  the  bow  and  stern  of  the  test  ship  to  a  tuned  bubble 
pulse  loading  (keel  shock  factor  1.19). 
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Appendix  -  A  sample  run  of  the  UBLOAD  program 


*  i|i  **  4  id  ♦  *  i|i  *  *  1(1  >ti  41  **  41  *  41 41 4141  *  41 1|>  *  *  41 41 41 41 1)1  *  *  i|i  id  41  i)i  i)i  41  *  1)1*  *  41  *  *  41 41  *  Id  41  *  *  iti  *  *  *  * 

Sample  run  ol  tha  program  UBLOAD.  Program  input  ia  in  lower 
case  lor  clarity.  *+*  indlcatea  comments. 
************************************************************ 


run  ubload 


ENTER  THE  6  CHARACTER  PREFIX  OF  THE  EQ.  BEAM  FILE 


ebcpl 

FILES  ARE  TITLED  ebcpl. TOT  AND  ebcpl. FRM 


UNITS  ARE  METERS, SECONDS, MEGAPASCALS 


% 


NO.  OF  HYDROSTATIC  SECTIONS  AVAILABLE  IS  20 


SCALE  FACTOR  IS  1.000000 


i  1  / 


COMPUTING  ADDED  MASS  VALUES. . . .HANG  ON 

***  Program  haa  located  the  required  prelx.tot  and  prelx.lrm 
+*+  files.  Bouyancy  and  fluid  added  matte  terms  are  being 
***  calculated. 

ENTER  THE  DEPTH  OF  DISCHARGE  (METERS) 


\Vl 

I  ,  • 


ENTER  PERPENDICULAR  STANDOFF  FROM  THE  SHIP  CENTERPLANE 


;*l 


***  A  value  ol  aero  Indicates  charge  is  directly  under  keel 
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ENTER  DISTANCE  FROM  FORWARD  PERPENDICULAR  TO  THE 
STANDOFF  PERPENDICULAR  -  NEGATIVE  IF  LOCATION 
FORWARD  OF  THE  BOW 

70 

ENTER:  1  -  FOR  SPECIFICATION  OF  KEEL  SHOCK  FACTOR 

2  -  FOR  SPECIFICATION  OF  CHARGE  MASS  (KG  TNT) 

2 


ENTER  CHARGE  MASS  IN  KG  OF  TNT 
167 


KEEL  SHOCK  FACTOR  IS  1 . 180070 


BUBBLE  MIGRATION  TO  BE  INCLUDED?  (Y-l.N-O) 

1 


FREE  SURFACE  EFFECTS  TO  BE  INCLUDED? <Y*1 ,N-0) 

1 

***  II  migration  and  Ira#  surface  effects  are  neglected,  the 
***  Tchebycheff  method  is  offered  as  an  option 

ENTER  DRAG  COEFFICIENT  CD  <0  FOR  DEFAULT) 

2.26 


***  A  default  value  of  2.26  is  used  in  the  program 
PULSE  NO.  1 

LSCALE  TSCALE  DELTA  ZINIT 
6.207467,  0.4162806,  1.610063,  6.477274 
ROOTS  ARE  0.1012467,  0.0330422 
BUBBLE  COEFFICIENT  K  IS  0.1703004 
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PULSE  NO.  2 

LSCALE  TSCALE  DELTA  ZINIT 
6.293110,  0.3540632,  1.850245,  6.423434 
ROOTS  ARE  0.2247307,  0.8706202 
BUBBLE  COEFFICIENT  K  IS  0.3260855 

K  VALUE  EXCEEDS  MAXIMUM  FOR  WHICH  ROOTS  ARE 
AVAILABLE  ON  PULSE  3.  NO.  OF  PULS.  •  LIMITED 
TO  2  FOR  THIS  EXPLOSION 

***  See  Section  2  of  the  report  for  a  discussion  of  the  k  value. 
***  Note  that  the  value  incraasea  from  pulae  to  pulae.  The  roota 
***  eorreapond  to  the  minimum  and  maximum  radii  predicted  by  the 
*++  deep  bubble  approximation  formula.  The  minimum  la  uaed  as  the 
***  initial  condition  for  the  integration.  The  acale  factora,  non- 
***  dimensional  depth  (DELTA) ,  and  non-dimensional  initial  head 
***  (ZINIT)  parameters  are  printed  for  reference  only. 

C0NTINUE7  (Y-l.N-O) 

1 

ENTER  8  CHARACTER  FILENAME  FOR  STORAGE  OF  BUBBLE  CHARACTERISTICS 
D  FOR  DEFAULT 
data . out 


ENTER  INTEGRATION  CUTOFF  TIME  IN  SECONDS 
ZERO  FOR  DEFAULT  VALUE 
3.0 


***  Typically  no  longer  than  2. 0-3.0  seconds  required. 

***  The  default  will  give  4*TSCALE,  or  almost  3  bubble  periods 

PULSE  TRANSITION  AT  TIME  STEP  113 

***  Transition  points  are  points  of  minimum  radii. 
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***  Energy  reduction  is  applied  at  those  points  if  required. 


BUBBLE  TIME  HISTORY  CHARACTERISTICS  WRITTEN  ON  FILE  data. out 


ENTER:  1  -  TO  CREATE  THE  VAST  .LOD  FILE 

2  -  TO  OUPUT  LOAD  TIME  HISTORY  DATA 

3  -  TO  OUTPUT  PRESSURE  TIME  HISTORY  DATA 

4  -  RESTART 
6  -  STOP 

NOTE:  ALL  FURTHER  OUTPUT  IS  WRITTEN  TO  FILE  ebepf .LOD 


WRITING  LOAD  FILE. . . 

ENTER:  1  -  TO  RESTART  ANALYSIS 

2  -  TO  OBTAIN  OTHER  DATA  FOR  THIS  RUN 

3  -  TO  STOP 

3 

CPU  tine  3.12  Elapsed  tine  1:61.86 


r 
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Figure  11:  Bubble  depth  time  histories  for  a  charge  of  227  Kg.  of  TNT  at  a  depth  of  45  m 


